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■ Abstract 

A formalism for the electromagnetic production of hypernuclei is developed where the cross sec- 
tion is written as a contraction between a leptonic tensor and a hadronic tensor. The hadronic 
tensor is written in a model-independent way by expanding it in terms of a set of five nuclear 
structure functions. These structure functions are calculated by assuming that the virtual photon 
interacts with only one bound nucleon. We use the most recent model for the elementary current 
operator which gives a good description of the experimental data for the corresponding elementary 
• process. The bound state wave functions for the bound nucleon and hyperon are calculated within 

vc : 

C\) . a relativistic mean-field model. We calculate the unpolarized triple differential cross section for 

the hypernuclear production process e + 12 C — > e + K + + as a function of the kaon scatter- 
ing angle. The nuclear structure functions are calculated within a particle-hole model. The cross 
section displays a characteristic form of being large for small values of the kaon scattering angle 
with a smooth fall-off to zero with increasing angle. The shape of the cross section is essentially 
determined by the nuclear structure functions. In addition, it is found that for the unpolarized 
triple differential cross section one structure function is negligible over the entire range of the kaon 
scattering angle. 
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I. INTRODUCTION 



The study of strange particles and hypernuclei remains an area of intense theoretical and 
experimental activity. Hypernuclei represent an exotic state of matter since they contain 
particles with quantum numbers such as strangeness which are not associated with ordinary 
nuclear matter. Whereas the nucleon-nucleon interaction is very well-known from elastic 
scattering data, our knowledge of the hyperon-nucleon, and hyperon-hyperon interactions 
is still relatively incomplete. A major goal of nuclear physics should also be a unified 
understanding of baryon-baryon interactions. Studies of hypernuclei are also important 
since they can give us insight into the role of strangeness in the stellar environment. A free 
A particle is unstable and will primarily decay via the weak interaction to a nucleon-pion 
system. However, a A in the nuclear medium will interact strongly with the other nucleons, 
hence forming a hypernucleus. The A is unaffected by the Pauli exclusion principle and 
can therefore occupy any one of the states already filled by the nucleons. In a one-boson- 
exchange picture the zero isospin of the A forbids the exchange of isovector mesons such 
as a pion or the rho meson with a nucleon, and therefore leads to a lack of strong tensor 
components in the AN interaction. Consequently the AN interaction is much weaker than 
the NN interaction, and the A in the nucleus does not lead to a major disruption of the 
shell structure jl|. Hypernuclei therefore provide remarkable experimental evidence for the 
shell-model of nuclear structure. 

A wealth of experimental data on hypernuclei have been accumulated by making use of 
hadronic probes [1]. These include the strangeness exchange reaction 

K~ (su ) + n (uud ) — > A (uds ) + tt~ (du ) (1) 

and the associate production process 

tt + (ud) + n(udd) — > A (uds) + K + (us). (2) 

An alternative production mechanism is through the use of real or virtual photons, i.e., 

>y + p(uud) — ► K + (us) +A(uds). (3) 

Hyperon production via the electromagnetic interaction requires the production of a strange 
quark/ ant i- quark pair. The large momentum transfer for associate production decreases 



2 



the sticking probability of the A, and consequently the probability for obtaining a bound 
hypernuclear system in the final state [2|. Since the production of the reaction particles are 
limited to the very forward angles, it is necessary to detect the electron and the kaon in 
coincidence jj}]. 

However, electron beams offer a number of distinct advantages. Indeed, with the advent 
of Jlab, our understanding of the role of the electromagnetic production process has greatly 
increased. The precision of the electron beam, as well as good spatial and energy resolution 
make up for the small (e, e', K + ) cross section relative to hadronic production mechanisms 
The (e, e', K + ) reaction converts a proton in a target nucleus, and populates proton-hole 
A-particle states. This reaction therefore produces neutron-rich A hypernuclei. Electro- 
production excites both natural and unnatural parity states with comparable strength (4]. 
Both the photon and the K + interacts relatively weakly with the nucleus and therefore the 
(e, e'K + ) reaction is not confined to the nuclear surface, hence hypernuclear states can be 
studied with a deeply-bound hyperon. In heavier nuclei the behavior of a A in nuclear matter 
may be studied. The transition operator has a spin-part, hence one can also probe spin-flip 
states. The electron beam can be polarized, whereas no polarized K~ (or 7r + ) beams exist 



The study of baryonic resonances is an important field in hadron phenomenology. Theo- 
retical work to determine the excitation spectrum of nucleons has been done mainly within 
the quark model framework. However, these models predict a much richer spectrum than 



hyperons |5j. In the electromagnetic production process, resonant baryon formation and 
kaon exchange play a primary role. The coupling of AT*'s and A's to meson-hyperon final 
states may be studied, and compared to SU(3) flavor symmetry predictions. In the case of 
electroproduction (q 2 7^ 0) two new features are introduced: (i) the longitudinal coupling 
of the photons in the initial state, and (ii) the electromagnetic and hadronic form factors of 
the exchanged particles. In the case of electroproduction, the cross sections for A and E° 
production is quite different. This is due to the isospin selectivity. In the K + A final states 
only the N* resonances are allowed, whereas for K + TP final states, A resonances may also 
contribute to hyperon formation. More information about the elementary process may be 
gleaned from electroproduction than from photoproduction, since the virtual photon mass 



|4j. 
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and polarization may be varied independently [6]. 

A number of experiments have been performed over the years to investigate strangeness ef- 
fects in nuclear physics. See for example Table I in Ref. Indeed, the experimental pursuit 
of the electromagnetic production of strangeness was given great impetus by the Jlab facility 
8]. Longitudinal and transverse cross sections were measured for the reaction 1 H( e, e'K + )A 
In experiment E89-009 the focus was on the production of hypernuclei [10|, lllf. The 
lypernucleus was produced via the reaction (e, e'K + ) using high-energy electron beams 
12|. Experiment E9 1-016 focussed on kaon electroproduction from deuterium [3, Q, as 



well as from 3 He and 4 He targets 



15. 



161 ] . The quasi- free electroproduction of unbound A, 
S° and X~ hyperons on carbon and aluminum targets was studied in Ref. |^]. Strangeness 
production off the proton and from nuclear targets has been investigated by the CLAS col- 
laboration. Cross section and recoil polarization data for the reactions 7 + p — > K + + A 
and 7 + p — > K + + S° for center-of-mass energies between 1.6 and 2.3 GeV 17|. 

The theoretical description of the elementary process is essential for studying hypernuclei 
formation. Obtaining results directly from QCD is a formidable task, and the standard 
approach is to use an effective field theory based on baryonic and mesonic degrees of freedom. 
The use of these so-called isobaric models has been pursued by a number of authors [5I. [l8 30] 
with a summary of theoretical work given in Ref. [311 ] . 

When making the transition from a theoretical description of the elementary process to 
that of hypernuclei production, a number of additional complications enter. These include 
(i) the description of the current operator in the nuclear medium, (ii) the nuclear structure 
model for the bound nucleons and hyperons, and (iii) the effect of nuclear distortion effects 
on the incoming and outgoing particles. The electromagnetic production of hypernuclei has 



)een investigated by a number of authors for photoproduction 



32- 



39j and electroproduction 



401 ] . The first complication is addressed by invoking the impulse approximation, i.e., the 



elementary current operator is assumed unchanged in the nuclear medium. The nuclear 
structure is described by non-relativistic shell model wave functions or by solving the Dirac 
equation with scalar and vector potentials to obtain bound state nucleon and hyperon wave 
functions. Finally, the nuclear distortion effects are treated within an optical potential 
formalism. In this work we adhere to the basic philosophy of these works, but also add a 
new feature, namely we write the triple differential cross section as a contraction of a leptonic 
tensor and a hadronic tensor. Following Refs. 



41 



431 ] we write the hadronic tensor in terms of 
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a set of five nuclear structure functions. Apart from the one-photon exchange approximation, 
and the additional assumption that the virtual photon interacts with only one bound nucleon, 
the formalism is still model-independent. In principle this allows one to study the cross 
section by doing a Rosenbluth-type separation. We do not explore this avenue in this 
work, but instead calculate the structure functions using a definite model for the process 
occuring at the hadronic vertex. Further, we calculate the bound state wave functions for 
the nucleons and hyperons using three different relativistic mean-field models. These are 
the linear Walecka model Q|, the successful NL3 parameter set J45], and the recently- 



introduced FSUGold parameter set 46|. Finally, we employ the most up-to-date form of the 
elementary current operator that gives a satisfactory description of the elementary process 
[ij]. In this work we use the model for the corresponding elementary process, together with 
the incorporation of nuclear structure effects, through an accurately-calibrated relativistic 
mean-field model. Thus, accurate binding energies and nucleon momentum distributions 
are employed. This method is based on the impulse approximation and provides a fully 
relativistic study both in the reactive content and the nuclear structure. In addition, it 
provides a simple way to study medium modifications of the produced mesons, as well as 
the resonances contributing to the elementary amplitude. 

This paper is organized as follows. In Sec. [Ill we discuss the kinematics, and show how 
the cross section may be written in terms of a contraction between the leptonic tensor and 
the hadronic tensor. Most of the technical details are deferred to Appendices |A] and [B] The 
model for the elementary amplitude as well as the nuclear structure models are discussed in 
Sees. Ill C 2l and lII C II Results for both the free process and the unpolarized triple differential 
cross section for hypernuclear production are given in Sec. II III with a summary in Sec. IIVI 

II. FORMALISM 

A. Cross section and kinematics 

Consider the electromagnetic production of hypernuclei 

e + A — > e + meson + aB (4) 
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FIG. 1: Lowest-order Feynman diagram for meson and hypernuclei electromagnetic production. 



which is shown schematically in Fig. [TJ If we confine ourselves to the extreme relativistic 
limit, then the incoming (outgoing) electrons may be specified by their four- momenta and 
helicity, i.e., (k,h) [(k',h')]. In the one-photon exchange approximation the reaction pro- 
ceeds via the exchange of a virtual photon with four-momentum g M = fc M — k ,fl = ((?o,q)- 
The outgoing meson is specified by its four- momentum p' x . The target and residual hyper- 
nucleus have four-momenta P and P', respectively. The differential cross section may be 
written in terms of these kinematical quantities and the transition matrix element M. as 

1 d 3 k' d 3 r>' d 3 ~P' 

da = ^3 9 V W V TT* (2?r)4 6{k + P ~ k ' - & ~ P ' } lM |2 (5) 

|vi - v 2 | [2tiY 2E p ' i {2'kY {2-k) a 

where |vi — Va | is the initial relative velocity. The transition amplitude A4 contains all 
the dynamical information about the reaction and will be studied in detail in Sees. Ill Bl and 
III CI In the laboratory frame the initial flux for massless electrons is equal to one. The 
spatial part of the four- dimensional Dirac delta function allows the integral over d 3 P' to be 
performed. This fixes the three-momentum of the residual hypernucleus to be 

P' = k-k'-pi = q-pi (6) 

where q = k — k' is the three-momentum transfer to the nucleus. In Appendix [X] it is shown 
that the triple differential cross section for the electromagnetic production of hypernuclei in 
the electron-nucleus laboratory frame is given by 

da 



dE k ,d (cos 9') dQ[ 



K\M\ 2 (7) 
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where K is a kinematic quantity that is fully determined by the energies and masses of the 
reaction particles, as well as the scattering angles of the ejectiles [see Eq. ( ]A14f) . 



B. Triple differential cross section in terms of leptonic and hadronic tensors 

Within the framework of the relativistic plane wave impulse approximation, the transition 
matrix element M for the electromagnetic production of hypernuclei may be defined as 



M = \U(k',h'h,U(k,h)] f£\ (pi; A *(P' )|>(g) |¥(P)> 



(8) 



with e 2 /47r = 1/137. In Eq. (jSJ) is the nuclear current operator, and U(k,h) is the 
plane wave Dirac spinor (defined in Eq. (1B2[) ) for the incident or ejectile electrons. \^{P) ) 
represents the many-body state for the incident nucleus, and \p\\ \E'(P / ) ) represents the final 
state consisting of the many-body residual hypernucleus state, and the outgoing meson. 
Using Eq. (JSJ) it follows that 

where we have introduced the leptonic tensor 



\M\ 2 = MM* = ( e - ) e^W" (9) 



£, u = [U{k',h') lfl U{k,h)] [U(k',h') lu U(k,h)]* (10) 
and the hadronic tensor 

W = (p'i, ${P') \ J^(q) \®{P) ) (p'i, i${P') \J u {q) \^{P) )*• (11) 

These two tensors are studied in detail in Appendix [B] In addition, it is shown in Appendix 
|B]that the unpolarized triple differential cross section for electromagnetic hypernuclei pro- 
duction in the electron-nucleus laboratory frame is given by 

± k (-X m wr 



,2 ) r Y S 

2 



dE k , d (cos 6' ) dn[ \q' 

= K (^) [W x (-3k-k' + 2f 1 (k,k'))+W 2 {-k-k'f 1 (P,P)+ 

\q 2 J E k E k i 

2h(k, P) fi(k', P)) + W 3 (-k ■ k' f x {p' 1: p\ ) + 2h{k,p' x ) f x {k',p\ )) + 
W, (2f 2 (P,p[))] (12) 

where K is a kinematic quantity that is fully determined by the energies and masses of 
the reaction particles, as well as the scattering angles of the ejectiles (see Eq. (IA14j) ). The 



functions /$ are defined in Eq. (1B37I) and (1B38j) . Apart from the one-photon exchange 
approximation and the additional assumption that the virtual photon interacts with only 
one bound nucleon, Eq. (IT2"|) is still model-independent. It shows that the unpolarized triple 
differential cross section may be determined from purely kinematical quantities, and a set 
of four nuclear structure functions, W\ to W4. In this sense Eq. (1121) is in line with the 



philosophy of Refs. 



42 



43 
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In the following section we present a calculation of these 



structure functions by evaluating the matrix element 



>' i;a *(P')I>(9)I^)> 



(13) 



in a model-dependent way. 



C. Model-dependent form of the hadronic tensor 

In the previous section a general formalism was developed for the electromagnetic pro- 
duction of hypernuclei. We now present a model-dependent evaluation of the structure 
functions Wt. 

The exact expression for the hadronic tensor is given in Eq. (1B6|) , and is defined in terms 
of the following matrix element (and its complex conjugate): 

J» = (p[; MP')\Hq)MP))- (14) 



To obtain a tractable form for this extremely complicated object, we rely on a number of 
approximations which are depicted schematically in Fig. |2j The principle assumption is that 
the virtual photon interacts with only one bound nucleon. This neglects two- and many- 
body components of the electromagnetic current operator. Secondly, it is assumed that 
the resulting meson and hyperons are produced from the interaction between the virtual 
photon and the nucleon to which it had coupled. This neglects two- and many-body rescat- 
tering processes. Additionally, nuclear distortion effects on the kaon are neglected. These 
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^ N meson ( p x ') 




meson (p x ') bound 
qs(P') \^ hyperonfp'. 




bound 
nucleon (p) 



FIG. 2: Graphical representation for the approximation employed at the hadronic vertex in order 
to obtain a tractable form for the matrix element given in Eq. (|13p . 



u ClA,CtN 



simplifying assumptions lead to the following expression for the hadronic matrix element 

d 3 p"d 3 P 5(q+ P - Pi - p") U« A (P") >(?)Z4 w (p) (15) 
d 3 pZ7 QA (q + p - pi) >(g) U aN (p) (16) 

(27T 2 iw) (p 2 sm6) Z7 aA (q+ p - pi) >(g) W Qjv (p)(17) 



where 



P = PmaxP 



= 7T0'' 



2tt <£". 



(18) 



(19) 



(20) 



In Eqs. (|T5l) - ( JTTjl the labels a a and o>n refer to the quantum numbers necessary to specify 
the bound state wave functions of the nucleon and hyperon, respectively. In Eq. (fl8~j) p max 
refers to the maximum momentum for which the momentum space wave function is still 
appreciable. More detail will be provided in Sec. Ill C II The operator J M refers to the current 
operator for the corresponding elementary process. The model for J M given in Eq. f[T7|) limits 
our results to simple particle-hole configurations of the produced hypernucleus. The model- 
dependent hadronic tensor is then given by 



^ model 



V (J" V 



(21) 



ajv,aA 



w 



(u- 1 ) u^WZ 



model 



(22) 



The model-dependent structure functions are then determined from the equation below (see 
also Eq. flB27j) ): 

W 2 
W 3 
\W 4 ) 

where the 4x4 matrix U is given by Eq. ([B29]). The two critical ingredients namely the 
bound state wave functions and the elementary current operator will be discussed in the 
following two sections. 



1. Nuclear structure 



In this work the bound state wave functions for the nucleon and the lambda are both 
determined from relativistic mean-field theory. For the nucleon bound state wave functions 



we use the linear Walecka model 



44j . the successful NL3 parameter set 45], and the recently- 



introduced FSUGold parameter set 46[ • For the hyperon state we employed the Lagrangian 
density of Ref. [48 1. 

For spherically symmetric nuclei the single-particle bound state wave function in position 
space is given by 



Z4(x) 



U EKm (yL) 



( 



\ 



9ek{t) 
r 

ifE K (r) 



y+«m(x) 



\ 



/ 



(23) 



where m is the magnetic quantum number, E is the binding energy, k the generalized angular 
momentum, and the spinor-spherical harmonics are defined as 

JWx) = ( l 2> m ~ S 'z> S 'z> ) Y l,rn-s> z {x) Xs' z - (24) 

4 = ±1/2 

The orbital angular momentum /, and the total angular momentum j, may be obtained as 
follows 



J 



1 

N - 2 



(25) 
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and 



I 



K, K > 

-1 — K, K < 0. 

The momentum space bound state wave function is defined by 

U EKm {v) = J d 3 pe~ ip x U EKm (x) 

^ 9ek(p) y Km (t>) ^ 



where 



and 



9ek(p) = / dr rg EK (r)ji(pr) 
Jo 



/e«(p) = sgn(/c) / dr rf EK (r) j 2 j-i(pr) 
Jo 

where ji(z) is the spherical Bessel function. 



(26) 



(27) 
(28) 



(29) 



(30) 



2. Elementary scattering operator 

The approximation depicted in Fig. |2] shows that the electromagnetic hypernuclei pro- 
duction process is essentially determined by the elementary process 

7 (virtual) + nucleon — > meson + hyperon. (31) 

In this work we invoke the impulse approximation and employ the model for J M as discussed 

in Ref. j^]. There are six different reaction channels which may be explored using this 
formalism, namely 

e + p — ► e + K + + A (32) 

e + n — ► e + K° + A (33) 

e + p — > e + K + + Y? (34) 

e + p — > e + K° + Z + (35) 

e + n — > e + K + + ^- (36) 

e + n — ► e + K° + S°. (37) 
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An application of the formalism in this paper will however, only be to the production of 
A hypernuclei. The electromagnetic current operator for each of these reaction channels is 
written as 



> = A(s,t,q 2 )Mt 



where 



MS 
Ml 

M5 

MS 



5 m 
7 rrv{ 



5 A* 
7 TOj 



1 



7 



7 5 (7^ - m 

±(p.q + p' 2 . q ) (2pT -<f)-\ (2pi • g - <? 2 ) (p" + j# ) 



7 5 ^3 = 7 5 • 97 A 

—^aXfivVl Q 7 9 



pT< 



In Eq. 



= 7 m s = 7 (pi 9 -Pi-qq^) 

- 7 5 < = 7 5 W-?V)- 

s and i are the usual Mandelstam variables defined as 
s = (q + p) 2 = (pi+p' 2 ) 2 
t = (q- Pi f = (p' 2 - pf 
U = (q ~ p 2 f = (p[ - pf ■ 



(38) 

(39) 

(40) 
(41) 
(42) 
(43) 
(44) 

(45) 
(46) 
(47) 



In Eqs. (145]) to (14?]) the four-momenta of the bound nucleon and bound hyperon are denoted 
by p and p' 2 , respectively. The four-momentum of the bound nucleon is defined as 



\M p — , p sin 9 cos 0, p sin 9 sin 0, p cos 9 



(48) 



where p, 9 and <ft refer to the integration variables defined in Eq. (fT7|) . The invariant 
amplitudes Ai are determined using an isobar model Feynman diagrams are written 
down for the s-, t- and li-channels for kaon electroproduction from the nucleon. The following 
resonances are included Su^SO), Pn(1710), Pi 3 (1720), L> 13 (1895), K*(892) and 7^(1270). 



III. RESULTS 



As a first application of the formalism developed in Sec. [TT] we calculate the unpolarized 
triple differential cross section for the hypernuclear production process 



e + 12 C 



e + K + + 12 B. 



(49) 
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This was the first hypernuclear spectroscopy experiment via electroproduction performed at 
Jlab |l2|. Before presenting cross section results we first need to investigate the two critical 
components that enter our formalism: (i) the elementary operator J M , and (ii) the nuclear 
structure input, namely the bound state wave functions for the bound nucleon and bound 
hyperon. 

The underlying elementary process for reaction (H9|) is 



e + p 



e + K + + A. 



(50) 



It can be shown that the triple differential cross section for the elementary hyperon produc- 



tion process can be written as 



25|: 



da 



K„ 



da 



V 



(51) 



dE k ,d{cosd')dQ[ ^ lep dtt[ 
where Ki ep is a kinematic factor that is purely described by the leptonic kinematics, and 
dav/dQ[ represents the differential cross section for kaon production from a virtual photon. 
This can be expanded in four terms which are related to the polarization of the virtual 
photon. In particular for longitudinally polarized virtual photons we have that dai/dVL^ 
may be written in terms of the /i = 3 and v = 3 components of the hadronic tensor 25], 
i.e., 



dai 



-2|p^ | y^g 2 
[{s-M*)ql 



E pi E p' 2 
2s 



W 33 . 



(52) 



Similarly the unpolarized transverse cross section is given by 

dQ[ 2{s-M 2 ) V > 
In Eqs. ( 1521) and (1531) the hadronic tensor for the elementary process is defined as 



(53) 



w"" = E J" (j v )* 

si,s' 2 



(54) 



E A A* Tr [Mf (fa + M) MJ (/ 2 + M Y ) ] 



(55) 
(56) 
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In Eqs. (I54j) - ( 1561) we have summed over the spin projections of the initial proton and the 
outgoing hyperon. For the elementary process the baryons may be represented by plane wave 
Dirac spinors U (p, s). The hadronic tensor was calculated in two independent ways: (i) from 
Eq. (|55[) by explicitly programming the Dirac spinors Z7(p, s) , the matrices Mf* (defined in 
Eqs. (139|) to dSD) and the current operator J M (defined in Eq. (155]) ). and (ii) from Eq. (|5"6"|) 
by explicitly performing the trace algebra over the Dirac matrices (i.e., the matrices Mf ). 
Both methods give identical results. This ensures that the current operator inserted in 
Eq. (iT7|) has been correctly implemented numerically. In Fig. [3] we show a graph of the 
longitudinal and unpolarized transverse cross sections as a function of Q 2 = —q 2 for the 
elementary process e+p — > e+K + +A with W = y/s = 1.84 GeV and the kaon scattering 
angle 8[ = 0°. The data are from Ref. The relatively good prediction of the scarce 
data by our model for the elementary current operator provides the confidence to embed 
this operator in the nuclear medium, for describing reactions on nuclei, and therefore obtain 
quantitative results for the triple differential cross section for hypernuclei production. 

In our formalism nuclear structure effects enter exclusively in terms of the momentum 
distribution of the bound nucleons and hyperons, and are calculated within a relativistic 
mean-field approximation. As was mentioned previously, there were three models that we 
considered, namely the linear Walecka model |44j, the successful NL3 parameter set [45], and 



the recently-introduced FSUGold parameter set [46] . In Fig. H] we show the results for the 
upper (g(r)) and lower (/(r)) radial wave functions in position space obtained from these 
three models for the ls 1//2 and lp 3 / 2 proton orbitals of 12 C. At the wave function level there 
is no real discernable difference between the different model predictions. The upper (g(r)) 
and lower (/(r)) proton wave functions (employing only the FSUGold model) together with 
the upper and lower lambda wave functions are displayed in Fig. [51 The momentum space 
wave functions are calculated from Eqs. (|29|) and (I30I) . The results are shown in Fig. [6] for 
the upper g(p) and lower (f(p)) momentum space wave functions. Once again there is very 
little difference between the models. We note that the wave functions are appreciable only 
for p m 0.6 GeV. This fixes the parameter p max referred to in Eqs. i fTTj) and ( {TBI . 

In Fig. [7] we show the radial momentum space wave functions for the proton and the 
lambda. For the proton wave function we only employed the FSUGold model. The binding 
energies for the proton and lambda for the different orbitals of 12 C and \ 2 B is shown in 
Table [TJ These binding energies are needed for the calculation of the magnitude of the 
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FIG. 3: Longitudinal (top figure) and unpolarized transverse (bottom figure) differential cross 
sections for the reaction e + p — > e + K + + A, as a function of Q 2 = —q 2 (four-momentum 
transfer). The kinematical quantities are W = y/s = 1.84 GeV and the kaon scattering angle 



0°. The data are from Ref. 



9]. 



three- momentum of the outgoing kaon. See Eq. (1A15|) . 

Next we display in Fig. [H] results for the unpolarized differential cross section (Eq. ( II 2|) ) 
as a function of the kaon scattering angle 8[, based on the following choice of kinematics: 
{E k , 9', E k >, (p'} = {3 GeV, 5°, 2 GeV, 0°}. For the hypernuclear production process given in 
Eq. (149]) there are four particle-hole transitions which may be studied within our simplified 
model-dependent form for (see Eqs. ( 1T5]) .( 1T6|) . ( 1T7|) and ( 123]) ). The cross section has the 
same behavior for all the possible transitions namely, large for small angles and a smooth 
fall-off to zero with increasing angle. The two upper cross sections (solid and dashed lines) 
correspond to the probe interacting with a lp 3 ^ 2 valence proton. Of these, the cross section 
is highest for a A in the Is 1 ' 2 shell. The two lower cross sections (dotted and long-dashed- 
short-dashed lines) are for a proton in the Is 1 / 2 shell. Again, the Is 1 / 2 A yields a higher 
cross section. 
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FIG. 4: Upper (g(r)) and lower (f(r)) radial wave functions in position space for the Is 1 / 2 and lp 3 / 2 



proton orbitals of 12 C. The solid and dashed lines represent linear Walecka model predictions [44j |, 
the long-dashed-short-dashed represents the NL3 calculation |45f] , and the dotted line represents 



the FSUGold model prediction 



46] 



It is also instructive to plot the model-dependent structure functions W\ to W4 as a 
function of the kaon scattering angle 8[. This is shown in Figs. [9] and [10] for each of the 
four possible transitions under consideration. It is clear from these figures that the shape 
of the cross section is determined by the structure functions. In Fig. [11] we plot the total 
cross section (indicated by the solid line), i.e., summed over all four possible transitions. 
The dashed line represents a similar calculation, but where the structure function W\ has 
been neglected. This graph suggests that W\ is negligible for 8[ > 16°, and makes a very 
small contribution for angles less than 16°. To quantify this we show in the bottom figure 
of Fig. [11] the quantity A defined as: 

(da ( do \ \ 

dE k td (cos 9') dn' x ~ \dE k , d (cos 6') dMj no m 



A 



da 



\ 



dE k ,d (cos 9') dQ[ 



x 100% 



(57) 
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FIG. 5: Upper (g(r)) and lower (/(r)) radial wave functions in position space for the ls 1//2 and 
lp 3 / 2 orbitals of 12 C and ^ 2 B. The solid line represents the ls 1//2 proton orbital of 12 C, the dashed 
line the Is 1 / 2 orbital of j^ 2 B, the long-dashed-short-dashed line represents the lp 3 / 2 orbital of 
12 C, and the dotted line the lp 3//2 orbital of ^ 2 B. For this figure we only employed the FSUGold 
model. 

as a function of 6[ . This extremely small difference illustrates that W\ is truly negligible over 
a wide angular range. The unpolarized triple differential cross section is therefore essentially 
just a function of three structure functions, namely W 2 , W 3 and W±. This could, in principle, 
allow a Rosenbluth-type analysis, similar to electron-proton scattering, to be performed for 
hypernuclei electromagnetic production to disentangle the various structure functions. 

Finally, we show in Fig. [12] the unpolarized triple differential cross section as a 
function of the kaon scattering angle 9[, for the kinematical set: {Ek,6',Ek',4>'} = 
{1.8 GeV, 1°, 0.5 GeV, 0°}. The calculation shown is for the total cross section, i.e., we 
have summed over all four possible transitions allowed within our simplified particle-hole 
model. As before the cross section is high for small values of the kaon scattering angle, with 
a smooth fall-off to zero as the angle increases. 
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FIG. 6: Upper (g(p)) and lower (f(p)) radial wave functions in momentum space for the ls 1//2 and 
lp 3 / 2 proton orbitals of 12 C. The solid and dashed lines represent linear Walecka model predic- 
tions [44J], the long-dashed-short-dashed represents the NL3 calculation [451 ]. and the dotted line 



represents the FSUGold model prediction 
IV. SUMMARY 



46] 



In this work a momentum space formalism was developed for the electromagnetic pro- 
duction of hypernuclei. The basic philosophy is to write the cross section as a contraction 
of a leptonic tensor and a hadronic tensor. The leptonic tensor is dependent on the helicity 
of the incoming and outgoing electron beams. We can therefore calculate fully, partially 
or unpolarized triple differential cross sections. The hadronic tensor is written in terms 
of five nuclear structure functions. The merit of writing the cross section in this way is 
that it could in principle allow a Rosenbluth-type separation to investigate the nature of 
the hadronic tensor. In this work we have not explored this avenue, but instead calculated 
the hadronic tensor based on the following model: it is assumed that the virtual photon 
interacts with only one bound nucleon in the nucleus, and that the elementary operator 
is left unchanged when it is embedded in the nuclear medium, i.e., we invoke the impulse 
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P (GeV) 

FIG. 7: Upper (g{p)) and lower (f(p)) radial wave functions in momentum space for the Is 1 / 2 and 
l p 3/2 orbitajg f and \ 2 B. The solid line represents the ls 1//2 orbital of 12 C, the dashed line 
the Is 1 / 2 orbital of ^ 2 B, the long-dashed-short-dashed line represents the lp 3 / 2 orbital of 12 C, and 
the dotted line the lp 3 ^ 2 orbital of \ 2 B. For this figure we only employed the FSUGold model. 

approximation. There are two critical ingredients when calculating hypernuclei cross sec- 
tions. Firstly, the model that is used to describe the elementary operator, and secondly, 
the model that is adopted to describe the nuclear structure. To calculate the elementary 
current operator we first expand it in terms of a set of six invariant amplitudes. These 
invariant amplitudes are calculated by writing down the Born diagrams, as well as the s- 
, t-, and -u-channel Feynman diagrams. The following nucleon and meson resonances are 
included namely, S n (1650), Pn(1710), Pi 3 (1720), £>i 3 (1895), K*(S92) and 7^(1270). In our 
model nuclear structure, which enter exclusively in terms of the momentum distribution of 
the bound nucleon and bound hyperon, are calculated within the framework of relativistic 
mean-field theory. We considered three models namely the linear Walecka model as well as 
the NL3 and FSUGold models. On the wave function level there is not a very big difference 
between the models and therefore we employed the FSUGold model for all our cross section 
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FIG. 8: Unpolarized triple differential cross section for the hypernucleus production process e + 
12 C — > e + K + + 12 ^B as a function of the kaon laboratory scattering angle 9[. The incident 
electron laboratory kinetic energy is = 3 GeV, the outgoing electron laboratory scattering 
angle is 9' = 5°, the outgoing electron laboratory kinetic energy is Ey = 2 GeV, and the angle 
between the leptonic and hadronic planes is <j)' = 0°. The solid line represents the calculation for 
a proton in the lp 3 / 2 orbital, and a A in the lp 3//2 orbital, the dashed line a proton in the lp 3 / 2 
orbital and a A in the Is 1 / 2 orbital, the long-dashed-short-dashed line a proton in the Is 1 /* orbital 
and the A in the lp 3 / 2 orbital, and the dotted line a proton in the Is 1 / 2 orbital and the A in the 
Is 1 / 2 orbital. 

TABLE I: Binding energies for the proton and lambda for the different orbitals of 12 C and ^ 2 B. 
The four numbers in square brackets for the proton orbitals are the predictions for the binding 
energy (in MeV) of the different models in the order QHDI, QHDII, NL3 and FSUGold. 



proton i2 C 



lambda \ 2 B 



Is 1 ' 2 [42.97, 38.99, 49.70, 39.19] 

lp 3 / 2 [16.17, 12.86, 16.04, 13.67] 



1 S V2 



[12.31] 
[1.11] 



20 



0.00 




16 24 32 
0/ (deg.) 



FIG. 9: Model-dependent structure functions W\ and W2 as function of the kaon scattering angle 
8[. The solid line represents the calculation for a proton in the \p 3 l 2 orbital, and a A in the lp 3 / 2 
orbital, the dashed line a proton in the \p 3 / 2 orbital and a A in the Is 1 / 2 orbital, the long-dashed- 
short-dashed line a proton in the ls 1 ^ orbital and the A in the \p 3 l 2 orbital, and the dotted line 
a proton in the Is 1 / 2 orbital and the A in the Is 1 / 2 orbital. 

calculations. As a first application of our formalism we calculated the unpolarized triple dif- 
ferential cross section for hypernuclear electroproduction from 12 C. A simple model for the 
transition matrix element was adopted which only included particle-hole transitions. The 
calculations indicate that the cross section is high for small values of the kaon scattering 
angle, and falls off smoothly to zero with increasing angle. The cross section for the four 
possible transitions has a specific structure. The cross section for transitions from the lp 3 / 2 
proton shell is larger than transitions from the Is 1 / 2 proton shell. In turn, for a specific 
proton shell, the cross section for the Is 1 / 2 lambda shell is higher. The individual model- 
dependent structure functions were also calculated. The results indicate that the shape of 
the structure functions is remarkably similar to the shape of the cross section. In addition, 
it is found that the W\ structure function is negligible over a wide angular range of the 
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0/ (deg.) 

FIG. 10: Model-dependent structure functions W3 and W4 as function of the kaon scattering angle 
9[. The solid line represents the calculation for a proton in the lp 3 / 2 orbital, and a A in the lp 3 ^ 2 
orbital, the dashed line a proton in the \p 3 / 2 orbital and a A in the Is 1 / 2 orbital, the long-dashed- 
short-dashed line a proton in the ls 1 ^ orbital and the A in the \p 3 l 2 orbital, and the dotted line 
a proton in the Is 1 / 2 orbital and the A in the Is 1 / 2 orbital. 

kaon scattering angle. This indicates that the unpolarized triple differential cross section 
is essentially just determined by three structure functions. This could, in principle, allow 
a Rosenbluth-type analysis to be performed for for hypernuclear electromagnetic produc- 
tion. There are many other questions which may addressed, such as the role of resonances 
in hypernuclear production compared to the free process or the role of spin observables 
as an additional tool to study how sensitive the hypernuclear cross section is to the ele- 
mentary operator. Our formalism also allows the study of possible medium effects on the 
resonances. Further improvements also need to be made to the calculation of the transition 
matrix element which in this work was based on a simple particle-hole model. However, we 
have established a model-independent form of the unpolarized and polarized cross sections 
in terms of nuclear structure functions. Improvements to the calculation of the transition 
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FIG. 11: The top graph shows triple differential cross section as a function of the kaon scattering 
angle 9[. The solid line represents the total triple differential cross section, i.e., where a sum over all 
four possible transitions has been performed, and the dashed line represents a similar calculation, 
but where the structure function W\ has been neglected. The bottom graph shows the quantity A 
defined in Eq. (J5j 



matrix element will therefore only impact on the hadronic tensor. 
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FIG. 12: Unpolarized triple differential cross section for the hypernucleus production process e + 
12 C — > e + K + + as a function of the kaon laboratory scattering angle 9[. The incident 
electron laboratory kinetic energy is E k = 1.8 GeV, the outgoing electron laboratory scattering 
angle is 9' = 1°, the outgoing electron laboratory kinetic energy is Ey = 0.5 GeV, and the 
angle between the leptonic and hadronic planes is 4>' = 0°. The solid line represents the total 
cross section, i.e., where we have summed over all four possible transitions within our simplified 
particle- hole model for the transition matrix element. 

Appendix A: Kinematics for electromagnetic hypernuclei production 



As is shown in Sec. Ill Al the differential cross section is given by 
1 d 3 k' d 3 p'i d 3 P 



da 



(2tt) 4 (5(A; + P- k'-p[ -P') \M\ 2 . (Al) 



|vi-v 2 | (2tt) 3 2E p; (27r) 3 (2tt) 3 

We can simplify Eq. OAip by employing the spatial part of the four- dimensional Dirac delta 
function to do the integral over the three-momentum of the residual hypernucleus which 
leads to 

5{E k + M A - E k , - E p{ - E P .) 



da 



cfk' (fp'i \M\ 2 . 



(A2) 



2(27r) 5 ^ 

In order to derive a triple differential cross section which may be compared to experiment, 
it is necessary to derive expressions that fully specify the four-vectors of the incoming and 
outgoing electrons, as well as the outgoing meson. The kinematical set-up for hypernuclei 
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FIG. 13: Leptonic and hadronic planes for hypernuclei electro-production. 

production is shown in Fig. [13j The direction of the virtual photon three-momentum defines 
the z-axis, i.e., 

z = A- (A3) 

|q| 

The unit vectors x and y define the leptonic plane in Fig. [13J 
The right-handed coordinate system is completed by defining 

y = z x x. (A4) 

In the leptonic plane the electron scattering angle is 8', and the direction of the incident 
electron with respect to the z is denoted by the angle a. In the hadronic plane the meson 
scattering angle is denoted by 8[. The hadronic plane makes an angle <ft with respect to the 
leptonic plane. The three-momentum of the incoming electron k (defined with respect to 
the coordinate system in Fig. [TBI is given by, for massless electrons, 

k = (insula) x + (E^ cos a) z. (A5) 

The energy transfer to the nucleus is given by 

g = E k -E v . (A6) 
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The three-momentum of the outgoing electron k' is given by 

k' = [E v sin(a + 6')} x + [E k , cos(a + 9')] z. (A7) 

Since the virtual photon direction defines the z-axis, it follows that q = (0,0, |q|), and 
therefore we can determine the angle a in Fig. [13] by demanding that q x = k x — k' x should 
be zero. This yields the following equation for the angle a: 

■ 2 El sm 2 6' 

sin a = — = =^ — — — -. (A8) 

E% + E% -2E k E k , cos9' 1 ; 

Geometric arguments show that the laboratory three-momentum of the outgoing meson is 
given by 



pi = JE P > - M 2 K [ (sin 9[ cos <p') x + (sin 9[ sin <p') y + (cos 9[) z ] . (A9) 



To proceed any further in specifying the kinematics, we need to determine the energy of the 
outgoing meson E p > . This is done by using the Dirac delta function in Eq. (1A2j) . The total 
energy of the residual hypernucleus is given by 



E P , = VP' +aM a = V (q - Pi f +aM a (A10) 
q - pi f + [M A - [M p - 4^) + [My - 4 y) ) ] 2 (AH) 

where E™ and E { P are the bound state energies for the nucleon and hyperon orbitals, 
respectively. These two quantities will be obtained from a relativistic mean-field model of 
the nuclear structure (see Sec. Ill C lj) . The quantity Ep> can be written as a function of the 
outgoing meson energy E p > . In addition, using the relations 

d 3 k' = 2itEl dE k ,d {cos9') , (A12) 



d 3 pi = yjEfy - M\ E p , dE p[ dn[ (A13) 

the triple differential cross section for the electromagnetic production of hypernuclei can be 
written as 



da E l J E k ~ M k 

V 1 \M\ 2 (A14) 



dE k , d (cos 6' ) (Kl'i 2 (2vr) 4 \f'(E. 
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where the energy of the outgoing meson E p > is obtained by finding the roots of the equation 
f(E p[ ) = E k + M A - E v - E p , - [q 2 + E 2 p , -M 2 K - 2|q| ^E%, - M 2 K cos^ + 

(m a -M p + M Y + E^ ] - e[ 



J B 



(A15) 



where 



|q| = [E 2 k + El - 2E k E k , cosfl'] 



(A16) 



We can see that the leptonic four-vectors, as well as the outgoing meson four- vector are fully 
determined if we specify the following kinematical quantities [E k , E k ', 9', fl[ = (9[, <p' ) ]. 



Appendix B: Leptonic and hadronic tensors 



The transition matrix element M. for the electromagnetic production of hypernuclei may 
be defined as 

m = \u(k',h') 7 ,u(k,h)] (£\ bk&iniHoiMP)) ( Bi ) 

withe 2 /47r = 1/137. In Eq. (IBlj) is the nuclear current operator. Since we are neglecting 
the electron mass with respect to the total energy, we employ the helicity representation 42 1 
of the plane wave Dirac spinor 



U{k, h) 



V2 



yh(j) h (k) ) 



(B2) 



with 



/ r 9' X -id,' ■ & \ 

Oh', i cos — - e v sm- 

9' 9' 

y 8 h/>1 e %<t> ' sin — + dw-x cos — 



(B3) 



where the unit vector k is specified by the polar and azimuthal angles, 9 and 0, respectively. 
This spinor is non-covariantly normalized to IPU = 1, and corresponds to the normalization 
adopted for the bound state spinors. See Sec. Ill CI It follows that |A4| 2 can be written as 
a contraction between the leptonic and hadronic tensors, i.e., 

,2\ 2 



\M\ 



(B4) 
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where 



and 



[ U(k', ti U(k, ti) } [ U(k', ti ) 7 „ U(k, ti) } ' 



>' i; A ^(P') |>(g) |*(P) ) (pi; A ^(P')\r(q) |*(P) >*. 



(B5) 



(B6) 



We will now study each one of these tensors in detail. The leptonic tensor may be written 
as a trace over Dirac matrices. This is done by using the identity 



U(k, ti) U(k', ti 



4:E k 



(B7) 



It follows from Eq. (1B7|) that the leptonic tensors will in general be dependent on the helicity 
of the incoming and outgoing electrons. We may therefore distinguish four different cases 
for the leptonic tensor. In the first case both the incident and outgoing electron beams are 
unpolarized. In this case we define 



f(0) 



Tr 



7u 



£ U(k,h)U(k,h)) lv ( U(k',ti)U(k',ti)\ 

h = ±l / W=±l / 



1 



(h^ k v + k^k u k ■ k g^u) ■ 



(B8) 
(B9) 



E k E k i 

Note that is completely symmetric in \x and v. In the second case the incident electron 
beam is polarized and the outgoing electron beam is unpolarized. In this case we define 



Tr 
1 



7m (U(k, h) U(k, ti) ) 7 , ( £ U(k',ti)U(k',ti ] 

k k ^ G^afS 



\h' = ±l 



(BIO) 



(Bll) 



2 ^ 2E k E k , 

where we adopt the convention e 0123 = +1 for the Levi-Civita tensor. Note that the lepton 
tensor now contains an anti-symmetric term due to the polarization of the incoming electron 
beam. If the incident beam is unpolarized and the outgoing beam is polarized then we define 



Tr 



7m Y, U(k,ti)U(k,ti)\ 7, {U(k',h')U(k',h')) 



\h' = ± 



^■va {k ,h ; k). 



(B12) 
(B13) 



In the final case the incoming and outgoing beams are polarized. Now 



»„„(*. h; k', ti ) = Tr [ 7(i (£/(k, h) U(k, h) ) 7 „ (U(V, ti) U(V, ti))} (B14) 



vv> _i_ 
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Next we turn our attention to the hadronic tensor W^. The definition in Eq. (1B6j) shows 
that this is an extremely complicated object since it contains exact many-body wave func- 
tions. However, can only be a function of the three independent four-momenta, namely 
g M , P M and p^ . Note that four-momentum conservation fixes 

P> = q + p-p'^ (B16) 

We can therefore expand W Mt/ in terms of a basis constructed from , g M , P M , p'^} . This 
is similar to the approach in Ref. [42| with the exception that a parity conserving electro- 
magnetic current forbids the presence of any terms linear in the Levi-Civita tensor. The 
expansion for then assumes the form 

W"" = W[ gT + W 2 q»q u + P^P V + W' A p'£j% + (P"q v + q^P v ) + 
^6 (Pl<f + (fVi ) + ^ (pTP u + P^pX ) + W % {P»q v - q»P u ) + 
W 9 (PU - <fvX ) + W w (p?P v - P% v ) • (B17) 

At this point contains ten independent nuclear structure functions. However, the 

imposition of electromagnetic current conservation, i.e., 

q^ u = q v W v = (B18) 

reduces the hadronic tensor to the following form 471 ] 



= W x + W 2 A»A V + W 3 B»B V + W 4 {A»B V + B"A V ) + 

W 5 {A^B U - B^A V ) (B19) 



where 



G nv = g »v _ y_y / B20 a 

T 

A^ 1 = P^ 1 — <f (B21) 
q l 

B» = p f - q*. (B22) 
T 

Note that 

A ■ q = B ■ q = 0. (B23) 
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The hadronic tensor consists of four terms which are symmetric with respect to fj, and v, 
and the last term which is anti-symmetric with respect to /i and u, hence 

4 



U 



[IV 



(B24) 



i = l 



where 



{<" ,u% u ,wf ,wf ,uf } = {G^A^S^A^ + S^A^-S^}. (B25) 

Since the contraction of a symmetric and an anti-symmetric tensor is zero, it follows imme- 
diately that 



IIV 



(B26) 



Since the basis {uf u } is not orthogonal, we can determine the structure functions Wi (i 
1, 2, 3, 4) by solving the following set of coupled linear equations 



Wjj = uw 



where 

(W u ) l = u^W, 
and the 4x4 matrix U is given by 



(B27) 



1,2,3,4 



(B28) 



U = 



( 3 


A 2 


B 2 


2A-B \ 


A 2 


A 4 


(A-Bf 


2A 2 (A ■ B) 


B 2 


(A ■ Bf 


B A 


2B 2 (A ■ B) 


\2A-B 


2A 2 (A ■ B) 


2B 2 (A ■ B) 


2A 2 B 2 + 2 (A - B) 2 J 



(B29) 



with 



W 



(W x \ 



(B30) 
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Using the general expansion of we can now work out its contraction with . As shown 
previously, the leptonic tensor can be written in four different forms, depending on whether 
the incident and/or outgoing electron beams are polarized or not. The four contractions are 

,2\ 2 



\M 



(0) 12 



- 2 ) e^r 

e 2 \ 2 1 



[W 1 (-3k-kf + 2f 1 (k,k?)) + W 2 (-k-k'f 1 (P,P)+ 

q J tik-bk' 

2A(fc, P) ft{k', P)) + W 3 (-k ■ k> ) + 2/ 1 (fc,pi) h{k',p\ )) + 



W, (2/ 2 (P,pi 



(B31) 



,2\ 2 



i-m (i) i 2 = ^ e^M')^ 



\ y / _ re re 



(B32) 



(B33) 



|.M (2) 



2\ 2 



2\ 2 



(B34) 



where 



l-M 



,2\ 2 



,2\ 2 



£^(k,h; k',h')W^ 



4 



(B35) 
(B36) 



/i(x,y) = x-y 



x ■ qy ■ q 

q 2 



h(x,y) = fi(k,x)fi( l kf,y) + f 1 (V,x)fi(k,y)-k-tff 1 (x,y). 



(B37) 
(B38) 



We can now substitute Eqs. flB3ll . flB33l) . flB34l) or flB36l) into Eq. fTATil) to obtain an 
unpolarized, partially polarized or fully polarized cross section. The resulting cross section 
will only depend on kinematical quantities (determined by the experimental set-up), and a 
set of nuclear structure functions. 
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